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We compute the perturbative expansion of the pressure of hot QCD to order .g 6 lng in the presence of 
finite quark chemical potentials. In this process we evaluate all two- and three-loop vacuum diagrams of the 
theory at arbitrary T and [i and then use these results to analytically verify the outcome of an old order g A 
calculation of Freedman and McLerran for the zero-temperature pressure. The results for the pressure and 
the different quark number susceptibilities at high T arc compared with recent lattice simulations showing 
excellent agreement especially for the chemical potential dependent part of the pressure. 
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1. Introduction 



If either the temperature or the density of strongly interacting matter is increased enough, it undergoes a 
phase transition from the hadronic phase into deconfined quark-gluon plasma (QGP) at an energy density 
of approximately 1 GeV/fm 3 . As one then approaches even higher energies, the value of the gauge coupling 
keeps decreasing making it possible to start using the machinery of perturbation theory in computing 
different observables. The problem of determining the perturbative expansion for the most fundamental 
thermodynamic quantity, the grand potential Q = —pV, has been under attack already for more than two 
decades. It is an especially hot topic today due to the fact that QGP is currently under experimental study 
in the ongoing heavy-ion experiments at RHIC. 

At vanishing chemical potentials, or zero net baryon density, the perturbative series for the pressure has 
recently been driven to the last fully perturbative order 2 , g 6 ln<7 ^JE], following the determination of the 
contributions of orders g 2 0, g 3 U, g^ng 0, g 4 and g 5 jZUHj. At zero temperature and large chemical 
potentials the expansion is known to 0(g 4 ) [S], and at high temperatures but finite chemical potentials to 
0(g 4 \ng) 0- The limit of large chemical potentials and small but non-zero temperatures is at present the 
least well known; there the only applicable result is of order g 2 In addition to these computations, there 
have been numerous attempts to determine the pressure using four-dimensional lattice simulations (see e.g. 
[n3IIIllI21IE3ini) and the HTL-approximation O [H [T7| . In the limit of a large number of flavors the 
pressure has furthermore recently been non-perturbatively determined both at /z = ^H] and \x ^ |19| . 
and these results have then been used to extract the perturbative expansion of the quantity at large n f even 
to order g 6 [T5] . 

The present paper provides a generalization of the order g 6 hag computation of Kajantie et al. to finite 
quark chemical potentials. Using the framework of dimensional reduction and evaluating all one-particle 
irreducible vacuum diagrams of the theory up to three-loop order, we will derive an analytic expression for 
the pressure valid at high temperatures and finite chemical potentials. Furthermore, we will determine the 
different quark number susceptibilities at fj, = 0, which together with the pressure are compared with recent 
lattice results of Gavai and Gupta showing impressive agreement in particular for the chemical potential 
dependent part of the pressure. The diagrammatic computations performed here will also be used to tackle 
the problem of determining the pressure at low temperatures. In particular we will verify the outcome of 
the well-known T — computation of Freedman and McLerran |2| and provide a simple analytic value for 
a poorly- known numerical coefficient appearing in the result. 

The paper is organized as follows. In section 2 the general notation is explained, and the necessary special 
functions are introduced. Section 3 provides then an introduction to dimensional reduction, and the results 
for the pressure at high T and finite /i are presented. These results are analyzed in detail in section 4, where 
we in particular investigate their agreement with lattice data. In section 5 we address the difficult problem 
of computing the pressure at T = and show how the result of [5] can be obtained from the computations 
performed in this paper. Section 6 is then devoted to addressing the important question of the compatibility 
of the two results obtained for the pressure at high T and T = 0, and conclusions are finally drawn in section 
7. We leave almost all computational details to the appendices. 

2. Setup and notation 

The theory we consider in this paper is the SU(AT C ) Yang-Mills theory coupled to nj flavors of massless 
fcrmions. It is described by the Lagrangian density 

£qcd = \F^F; V + ^, (2.1) 

2 At order g 6 one runs into infrared problems that can only be solved by non-perturbative methods 1201 . 
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where, as usual, 



= d^-d^+gf^A^, (2.2) 
Dp = d^-igA^ = d^-igA;T a , (2.3) 

and the symbols T a denote the generators of the fundamental representation of the gauge group. All quark 
fields have been combined into a multi-component spinor ip, and since flavor is a conserved quantum number 
of the theory, we may assign an independent chemical potential \if for each /. As is customary in finite 
temperature computations, we will work in Euclidean metric. 

In finite temperature field theory the partition function is represented by a functional integral of the 
exponential of the Euclidean action, where the usual time integral has been replaced by one over the 
compact imaginary time r (= x a in Euclidean metric) ranging from to l/T 

Z(T,n) = e- n ' T = yp^cxpj- dr J ^-^(C- MO}- ( 2 - 4 ) 

The perturbative evaluation of this integral leads to the computation of vacuum Feynman diagrams with 
Feynman rules closely analogous to the zero-temperature ones. The most important modification is the 
replacement of the p loop integrals by discrete sums over the so-called Matsubara frequencies 

3^ os = 2nnT, (2.5) 
pf° rm = (2n+l)7rT-z M , (2.6) 

where n is an integer. In gauge field theories such as QCD the gauge invariance creates an additional 
problem, as one needs to restrict the degrees of freedom contributing to the functional integral to the physical 
ones. In the present paper this is implemented by working in the covariant Feynman gauge throughout the 
computations. 

We end the section by introducing some new notation. The chemical potentials will henceforth usually 
appear in the dimensionless combinations 

ft = m/(2ttT), (2.7) 
z = 1/2 -ip, (2.8) 

and in the context of computing the zero temperature partition function the following abbreviation will be 
used 

£"/ = (2-9) 
/ 

In sums over a single flavor index the subscript / in [if will usually be suppressed. 

The momentum integration measure and the notation used for sum-integrals from here onwards are 



f f ddp A - 2e pA 2 V f d d p 

L J (2n) d " I 4^ J J (2nr ' 



(2.10) 



p/{p} 

Po 



2 / 



t yi, (2.ii) 



where A is the MS scale, and po and {po} denote the bosonic and fcrmionic Matsubara frequencies, respec- 
tively. We have introduced the unconventional notation A for the MS scale in order to avoid confusion with 
the chemical potentials. 
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Wc define some familiar group theory factors by 

C A S cd = f abc f abd = N c 5 cd , (2.12) 

Nj_l 
2N B 



C F b l3 ee (T a T%- = ( 2 - 13 ) 



T F 5 ab ee Tr T a T b = ^ S ab (2.14) 
and denote an additional, slightly less well-known one by 

N 2 — 4 

DS cd = d abc d abd = — S cd . (2.15) 

N c 

The dimensions of the adjoint and fermionic representations of the gauge group are naturally 

d A ee 6 aa = N^-l, (2.16) 

d F ee S u = d A T F /C F = N c n f . (2.17) 

For some frequently occurring combinations of special functions we will apply the following abbreviations 

C(x,y) ee 8 x C{x,y), (2.18) 

H(n,w) ee C(-n,w) + (-l) n+1 ('{-n,w*), (2.19) 

K(w) ee *(«;) + *(«;*), (2.20) 

where n is assumed to be a non-negative integer and w a general complex number. Here £ denotes the 
Riemann zeta function, and \P is the digamma function 

*H - (2.21) 
These functions are analyzed in some detail in appendix D. 



3. The pressure at large T/fi 
3.1. Dimensional reduction 

In order to compute the partition function of QCD, we need a systematic way of taking into account the 
contributions of the different momentum scales, as conventional perturbation theory fails already at three- 
loop order. At high temperatures a physically intuitive solution is offered by dimensional reduction, which 
is based on the observation that as the temperature is increased, all degrees of freedom except for the 
ones associated with the zero Matsubara modes of bosonic fields get large effective masses proportional to 
T. They can thus be integrated out leaving us with a three-dimensional effective theory describing the 
soft scales, where only the bosonic zero modes remain intact. Details of the procedure can be found from 

The inclusion of the chemical potentials makes the problem of determining the pressure even more com- 
plex, as we now need to take into account the effects of these new scales in addition to the usual thermal 
ones proportional to T. Assuming the magnitude of the chemical potentials to be negligible in comparison 
with 2ttT we may, however, certainly continue using dimensional reduction as a framework |23| . This means 
that the expression for the pressure may be separated into three parts 

PQCB = Pe+Pm+Pg, (3.1) 
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corresponding to the contributions of the momentum scales 2ttT, gT and g 2 T, respectively By 
definition 



p E (T, n) = PQ cB(T,n)-ylnJvA«VA a ex V {-S E }, 
where Se is the action of a three-dimensional effective theory with the Lagrangian density |23| 
C E = \^F?j + Tr[D l ,A } 2 +m 2 E TrA 2 + ^ ^M/Tr^ + <5£ E , 



(3.2) 



(3.3) 



and where the traces are now taken only over the color indices. Similarly, pm and pg are defined by 



Pm(T, m) 



Pqcd(T, h)-Pe(T,ix) 



T 
V 



In /" £U?exp{-S M } 



PQCT)(T, n) - p E (T, fi) - p G (T) , 



C M = -TrF^ 



<5£ 



M • 



(3.4) 
(3.5) 



The gauge coupling constants of the two effective theories, <?e and gM, appear in the covariant derivatives 
above, and operators contributing to the partition function starting at 0(g 6 ) or higher have been assembled 
to the terms SCe and 5 Cm . The question, at which values of the chemical potentials we may trust results 
obtained using dimensional reduction, is examined quantitatively in |23| and is also briefly discussed in the 
section 6 of this paper. 

At leading order the different parts contribute to the pressure as pe ~ 5°, Pm ~ g 3 and po ~ <? 6 lng. The 
first of these functions can be obtained by computing the strict perturbation expansion of the pressure in 
the full theory, i.e. by evaluating all 1PI vacuum diagrams of four-dimensional QCD without applying any 
form of resummation. The two other ones are then available by evaluating the perturbative expansions of 
the partition functions of the effective theories the parameters of which must, however, be determined 
through the full theory |23U24| . Following the notation of ^ and using results from 0123 1 these statements 
can be summarized by writing 
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Figure 1: The two- and three-loop fermionic diagrams of the full theory contributing to the values of «E2 
and «E3- The solid, wiggly and dotted lines stand respectively for the quark, gluon and ghost propagators. 
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(3.10) 
(3.11) 

(3.12) 



where g is the renormalized gauge coupling of quantum chromodynamics and the matching coefficients a are 
left to be determined. Apart from modifying the values of these coefficients, the effects of finite fi can only 
be seen in the appearance of the last term in Eq. |3.9p . In particular, one should notice that j?g depends 
only on T at the order considered here. 

With the exception of a?E3 and cues, the results for the matching coefficients can be immediately extracted 
from [Tll3*l l23U25| . To get «E5 we furthermore merely need to evaluate the one- loop gluon polarization tensor 
to 0(e) in the limit of vanishing external momenta, which is a simple computation. The calculation of «E3 
is, on the other hand, already a considerably more laborious task, as it involves computing all three-loop 
vacuum diagrams of the theory. They have so far only been evaluated at vanishing chemical potentials [H] 
and temperatures 9 , and the generalization of these calculations to finite fi and T is the topic of appendices 
A and B. The relevant fermionic two- and three-loop diagrams are depicted in Fig. 1, and the results for 
the matching coefficients can be found from below. 

In order to write the perturbation theory result for the pressure in the familiar form of a power series in 
the coupling constant, we simply need to add together Eqs. (|3.6|) . I|3.9|l and l|3.12|l and expand the result in 
g. Up to order g 6 hig the outcome reads 
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where the pole of a^3 exactly cancels the 1/e term appearing in the order g A contribution. As there is an 
unknown 0(g e ) contribution missing from the result, there is an ambiguity in choosing the coefficients inside 
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the logarithms of the order g 6 lng terms 3 . The current choice is, however, well-motivated from the effective 
theory point of view; the logarithms now appear in the same form they emerged from the three-dimensional 
Feynman diagrams. 

3.2. The matching coefficients 

Given in terms of the special functions and group theory factors defined in the previous section, the results 
for the matching coefficients a read 
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3 Varying the coefficients inside the logarithms amounts to varying the magnitude of the undetermined order g 6 term. 
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<*» = ^E{^( 221n ^+ 5 )+^( 2 ( 7 + 132 M 2 ) ln £^+ 9 + 132 ^ 



8 7 + 4H(z) 



- 18 C F T F (l + 12^ 2 ) -4T F (l+ 12/2 2 ) ^2 ln^ - 1 - N(z) j j , (3.19) 
a E7 = ^E{^( 22 ^+l)-4^( 2 ln^-H(,))}. (3.20) 
Combined with Eq. (|3.13() . this is the main result of the paper. 



4. Lattice tests 

As there unfortunately is no experimental data available for the pressure in the QGP phase, the results 
derived for this quantity in the previous section can only be compared to other analytic computations 4 or 
to lattice simulations. In particular, it is very interesting to investigate, to what accuracy we can reproduce 
the results of the various lattice studies that have been performed for the pressure and the different quark 
number susceptibilities. In this section these comparisons will be made, and it will furthermore be studied, 
how rapidly the perturbative series for the different quantities converge. 

4.1. The pressure 

It has been observed that at vanishing chemical potentials the order <? 6 ln g perturbative result for the pressure 
is well compatible with four-dimensional lattice simulations, but that the eventual determination of the yet 
unknown g 6 term may still change the situation dramatically pQ. The perturbative result varies largely from 
order to order and even at 0{g ) its behavior as a function of temperature still bears no resemblance to the 
lattice predictions, even though the order g 2 result gives a relatively good estimate for the quantity. For 
the fj, dependent part of the pressure we will see a significant improvement in these convergence properties. 
Extracting the quantity 

AP(T,/i) = p QCD (T, M )-p QCD (T,0) (4.1) 

from Eq. (|3.13() , the /i dependence of the pressure can be directly compared with recent lattice studies ^1] , 
where AP has been computed in quenched QCD assuming two light flavors of quarks, u and d, at equal 
chemical potentials. In Fig. 2 (left) this lattice data is plotted alongside with the perturbative result, which 
has been obtained by setting all explicit factors of rif to zero in order to match the quenched approximation. 
One observes that already at temperatures T ss 2T C the perturbative results lie well within the error bars 
of the lattice datapoints and that the differences between subsequent perturbative orders are very small. 

As we can see from Fig. 2 (right), the picture is qualitatively similar to the fi — case in the sense that 
the leading correction to the free theory result already gives a good estimate for the quantity in question. 
The next perturbative orders then make the situation worse until at 0(g 5 ) one again starts approaching 
the lattice results. The main difference between the two cases is simply that AP is a much more strongly 
perturbative quantity: for it even the free theory result falls within 10% of the lattice data and the relative 
magnitudes of the perturbative corrections are considerably smaller than for the /i = pressure. 

4.2. Quark number susceptibilities 

Apart from analyzing AP directly, there are other, more effective ways to investigate the chemical potential 
dependence of the results. To make full use of the large amount of lattice data existing at fj, = (see e.g. 

4 A recent large-n^ computation by Ipp and Rebhan provides an accurate numerical check for some of the results of this 
paper. For details, see section 3 of 1191 . 
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Figure 2: The perturbative and lattice [T$\ results for AP plotted as functions of T/T c . On the left (i has 
been given different values, and on the right the different perturbative orders for the (i = 0.44T C case are 
shown normalized to the free theory result. The result T c /Aj^g|„ /= o = 1-15 |26j and the /i = convention 
A = 6.742T have been applied here. 



^^1231221 EDI) we ma y use Eq. I|3.13|l and the results of appendix D to compute the different quark number 
susceptibilities 

d n p 

Xijk --- d fad fijd (Ak- 

in this limit. The linear (i.e. second order) susceptibilities have already been considered both in the frame- 
work of ordinary perturbation theory [21] and in the HTL approximation |25l 1311 l3~2"] with the result that 
only the diagonal ones are accurately predicted by perturbation theory. For the linear non-diagonal suscep- 
tibility the perturbative results were found to be several orders of magnitude larger than the lattice ones 
25 . One should, however, note that for the non-diagonal susceptibilities there is considerable disagreement 
between the different lattice approaches (see the discussion in |30j ) . 

A similar behavior can be observed when studying the nonlinear quark number susceptibilities. In Fig. 
3 we have plotted the susceptibilities Xuuuu and Xuudd at fif — 5 , and comparing with the quenched QCD 
results of ^3] we again see that the diagonal quantity is satisfactorily produced by perturbation theory but 
that the prediction for the non-diagonal one is too large by more than a factor of 1000. This apparent 
disagreement is, however, not unexpected, as even the different lattice results for the non-diagonal suscep- 
tibilities differ from each other. Furthermore, the perturbative expansions for the nonlinear susceptibilities 
start only at relatively high orders (0(g 6 \ng) for Xud, C(<? 3 ) for Xuudd), and it is therefore entirely possible 
that large cancellations will occur as one drives perturbation theory even further. The situation is com- 
pletely different in the case of the the diagonal susceptibilities, for which the free theory result already gives 
the correct order of magnitude of the results. E.g. for Xuuuu one obtains from Eq. 13.14fl 

Xuuuu = pj = ^+0(g 2 ) ~ 0.61 + O(.9 2 ), (4-3) 

U(l u 7T 

which is in good agreement with the lattice data (see Fig. 3. a). 

5 In this context the rif = limit of the perturbative result is understood to be taken only after the necessary differentiations 
of the pressure with respect to fif have been carried out. The pressure at rij = is defined in an analogous manner. 
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Figure 3: Perturbative results for Xuuuu and Xuudd at n/ = plotted with lattice results from |14|. The 
curves corresponding to different values of A show the expected effect of the yet undetermined g 6 term of 
the perturbative expansion (for details, see |27]). Again T c /A^jg|„ /= o = 1-15 and A = 6.742T |24| . 



For AP, and hence for the susceptibilities, the order g 6 term in the perturbative series contains no 
non-perturbative contributions and is therefore in principle straightforwardly obtainable, unlike the corre- 
sponding term in the expansion of the \x = pressure pQ. This computation may already be enough to 
improve the perturbative predictions for the non-diagonal susceptibilities significantly, but will have practi- 
cally no effect on the already good convergence properties of the chemical potential dependent part of the 
pressure. As has been pointed out in ^3] > the effects of the nonlinear susceptibilities on AP are negligible for 
small values of the chemical potentials, and the quantity is almost solely determined by the linear diagonal 
susceptibilities. 



5. The pressure at T = 

The zero-temperature pressure of QCD was first computed to 0(g 4 ) a long time ago [5] in a lengthy 
calculation involving numerical integrations. Using the analytic results for three-loop diagrams derived 
in the present paper this result can be straightforwardly analytically reproduced. The computation is 
divided into two distinct parts: the results for the graphs that remain infrared convergent at T — may be 
immediately continued to this limit, but in addition an infinite set of IR divergent ring diagrams must be 
summed over explicitly. Analogously to the use of the three-dimensional effective theories in section 3, this 
resummation is necessary to ensure that the contributions of all momentum scales are properly accounted 
for. Only the results of the computation are given below, while the details are left to appendix E. 

5.1. IR convergent diagrams 

At T = QCD pressure gets contributions only from the fermionic graphs, i.e. from the diagrams of Fig. 
1. Aside from the IR divergent graph I e we obtain using the results of appendices A - D 

Pi = {otm + I a + h + I c + Id + If + I g + h)\ T=0 
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Figure 4: a) The fcrmionic part of the one-loop gluon polarization tensor divided into its vacuum and matter 
parts. 

b) The diagram I' e contributing to the zero temperature pressure. 

c) The generic form of the ring diagrams contributing to p^. 
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(5.1) 



where the renormalization of the gauge coupling has been taken into account. 



5.2. Ring diagrams 

In order to obtain the correct expression for the zero-temperature pressure up to order g A we need to add 
to Eq. (|5.1(l the contributions of all ring diagrams of the type I e . Individually these graphs are infrared 
divergent but when summed together they give a finite contribution to the pressure starting at O(gHng). 
Separating the fcrmionic part of the one-loop gluon polarization tensor into its vacuum (T = n = 0) and 
matter (vacuum subtracted) parts as in Fig. 4. a. we observe that to order g 4 only the diagram I' e of Fig. 
4.b. and the ring sum of Fig. 4.c. need to be computed. The reason for this is that starting at four-loop 
order the diagrams with at least one vacuum insertion only contribute at 0(g 6 hig) or higher, and the 
corresponding three-loop diagram with two vacuum insertions naturally vanishes at T = 0. 

A straightforward computation performed in appendix E shows that the diagram I' e gives the following 
contribution to the pressure 

4 



P2 

■4/7"" 

/ 

The summation of the ring diagrams was, on the other hand, first performed in [5] and is reproduced in 
appendix E following in most parts the treatment of the original paper. The result of this computation 
reads 

n = -£<,^{4,ni-f + M hSa _ M+ , + £ 



3 ^-(f 2 ) 2 (It 2 ) 2 



where we have defined 



4 M/Mg (/4 + /flln ^ + J*? - 04 - Mg )ln ^ | ■ (5.4) 



10 



The constant S possesses the integral representation 
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~ -0.85638320932694280684831023291594035884727909711135760899309086726768550829, (5.5) 

which we have not been able to evaluate in closed form. Even all attempts of expressing its numerical value 
in terms of the most common natural constants using the PSLQ algorithm j^S] have been unsuccessful. 
Finding the correct basis of constants for 6 seems to be a very non-trivial task. 



5.3. The result at T = 

Eqs. I|5.1|) i H5.2f) and (|5.3|l together verify the well-known result of |§j 



Pqcd{T = 0) 



Pi+p 2 +P3 + 0(g 6 \ng) 



I (UN c -2n f )hx- 



!6, „ 
— In 2 

3 



17 1 



1 



+ TTT +^ (415 - 264 In 2) A^ c - - — -In 2 n 



4 Nr. 36 

.4 



11 



(i 



47T, 



4 In- 



47T 



22 

y 



16 



2tt 2 



s hi2(l-m2)+£ + , 




In particular, we have here obtained an analytic value for the coefficient of the term, which was 
previously known only numerically with considerable error bars. When comparing Eq. I|5.6[) with the result 
of jH], one should notice that there the authors work in the momentum subtraction scheme, in which the 
gauge coupling constant is related to the one of the MS scheme through the equation 



4tt 2 
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4tt 2 



(5.7) 



6. Compatibility of the results at large /i/T 

In the sections 3 and 5 of the present paper we have derived perturbative results for the QCD pressure in 
the two limits of high temperature and small chemical potentials and T = and large \x. Determining the 
exact region of applicability for the first one is a nontrivial task, as it has not been studied analytically, 
how the appearance of the additional scales /i/ affects the validity of dimensional reduction. In |23| it has, 
however, been estimated based on numerical results for correlation lengths that the method applies as long 
as \if < 4T for all flavors. This seems physically very reasonable, since one certainly expects the framework 
of dimensional reduction to be unaltered, if the values of the chemical potentials are much smaller that the 
thermal scale 2-7rT. Perhaps not surprisingly, roughly the same result was obtained in |14) for the quantity 

y = — t 6 - 1 ) 

J V ± Xuuuu 

describing the highest value of chemical potentials, for which the linear susceptibilities accurately produce 
the /^-dependent part of the pressure, AP. In Fig. 5 (right) the perturbation theory and lattice results for 
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T 
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Figure 5: On the left a sketch of the \i-T plane, where F-McL stands for the T = result, Eq. I|5.6|) . Next 
to it the perturbative result for the quantity /i* plotted for n/ = together with lattice data from |14| . 
Once again T c /k m \ nf = = 1.15 26 and A = 6.742T gj. 

/x* are plotted as functions of temperature showing reasonable agreement. The overall scale of the results 
is again given by the free theory expression, which now reads 

^ = V2n + 0(g 2 ) ~ 4.4 + 0( 5 2 ). (6.2) 

There remains a region on the (i-T plane between the lines T = and T = /i/4, where the perturbative 
expansion of the pressure is only available to 0(g 2 ). To obtain an order g 4 result valid throughout the 
deconfined phase, one would have to perform an explicit summation of all the bosonic and fermionic ring 
diagrams at an arbitrary temperature, as adding a mass term for the zero mode of A Q to the free Lagrangian 
of the theory would not lead to the expected result in the limit T — > 0. Even though this procedure is enough 
to produce the correct 0(g 4 ) result for the pressure at high T it does not work at low temperatures 
due to the nontrivial structure of the gluon polarization tensor at a vanishing temperature and external 
momentum (see appendix E). The separation of the zeroth Matsubara mode of the Ao field is furthermore 
clearly inconsistent, if fj, ^> ttT. 

The interesting limit of small but non-zero temperatures can in any case be formally taken also in Eq. 
(|3.13|l , even though it is already beforehand understood that an unphysical logarithmic divergence of the 
type ln(T/A) will appear there. As T approaches zero, it is natural to investigate, how the difference of Eqs. 
(|3.13|) and l|5.6|l behaves as a function of the chemical potentials, as it gives the magnitude of the terms 
that have been neglected in deriving Eq. 1)3. 13[) but are necessary to obtain the correct T = pressure. 
This quantity has been plotted in Fig. 6 (right) for different temperatures alongside with the corresponding- 
curves for the pressure. In this figure we assume two flavors of massless quarks at equal chemical potentials, 
and the scale parameter has somewhat arbitrarily been chosen to be 

A = 2 V T2+ 2^ (6 ' 3) 

in analogy with the free theory pressure, Eq. (|3.14|) . 

The result shown in Fig. 6 is remarkable. It seems that as we approach the zero-temperature limit, the 
curves corresponding to Eq. (|3.13|) smoothly approach the one describing Eq. I|5.6|l . until only at very low 
temperatures T <C T c the logarithmic divergences start increasing the gap. This observation suggests that 
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Figure 6: On the left the behavior of Eq. (|3.13|1 as a function of the chemical potentials is portrayed for 
different temperatures alongside with the T = result, Eq. Ij5.6|l . On the right the difference of these 
results is shown on a logarithmic scale appropriately normalized. The result T c /Aug| n/=2 = 0.49 has 
been applied here. 



the magnitudes of the terms not present in Eq. (|3.13|) are small and one is able to use this result throughout 
the deconfmed phase with the exception of a narrow strip near the T = line. We in particular notice that 
a large value of fi/T does not itself appear to spoil the applicability of Eq. I|3.13[l an d that the restriction 
A* ^ 4T may therefore perhaps be loosened. 

Despite all the optimism, we must be very careful in interpreting Fig. 6: only one special configuration 
of chemical potentials has been analyzed so far, and the reasoning presented above is merely of qualitative 
nature. It is furthermore clear that dimensional reduction cannot be reliably applied in the limit of small 
temperatures. The good compatibility of Eqs. 1)3. 13[) and l|5.6|) is most likely simply a consequence of the 
fact that in both of them the numerically dominant part comes from the strict perturbation expansion of 
the pressure, Eqs. (|3.14() - l|3.16[) . Another aspect to keep in mind is that perturbative computations such 
as the one presented in this paper can never produce the rich phase structure of the 'condensed matter 
QCD' found in this region of the fi-T plane. Thus the applicability of the present results is in any case 
very limited there. 



7. Conclusions 

In this paper we have improved the perturbative expansion of the pressure of hot and dense QCD by three 
orders, as the g 4 , g 5 and <7 6 ln g terms in the series have been determined. The crucial step in the computation 
was the analytic evaluation of all three-loop 1PI vacuum diagrams of the theory at arbitrary T and [A, which 
were also used to derive an order g 4 result for the zero-temperature pressure. Finally, it was argued based 
on a qualitative analysis of our results that the perturbative expansion of the pressure is now converging 
relatively well on almost the whole /i-T plane. 

There is, however, a large amount of work left to be done. At high temperatures and small chemical 
potentials one clearly needs to determine the next 0(g e ) term in the perturbative expansion, as this order 
contains the first non-perturbative contributions to the pressure and furthermore has a potentially very 
significant impact on the behavior of the result. Another challenge can be found in improving the present 
embarrassing record in the limit of large chemical potentials and small but non-zero temperatures; a gen- 
eralization of the order g A result at T = to non-zero temperatures would certainly be welcome, even if its 
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numerical effect on the present results turned out to be small. 
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Appendix A. Vacuum diagrams 



In order to complete the calculations of section 3, we still need to evaluate the two- and three-loop 1PI 
vacuum diagrams of QCD. Since the purely bosonic graphs are unaffected by the finiteness of the chemical 
potentials and have been computed already previously 0, one may restrict the treatment here to the 
fermionic diagrams of Fig. 1. The calculations will be performed keeping both the temperature and 
the chemical potentials arbitrary which together with the relations listed in appendix D enables one to 
immediately continue the results to the limits of small T and small (i. 
Let us start by defining a set of 'master' sum-integrals 

xr s (A - 2) 



T = % (A.3) 

J{pq} P 2 Q 2 (P- Q) 

T = Y — (A 4) 

T {PQ} P 2 Q 2 (P - Q) 4 ' 



Mm ' n ~ f {PQR} P 2 QHRr i ((P-Qf) n {P-R) 2 {Q-R?' {A ' 5) 

Mm ' n " f {PQ }R PiQ* (i? 2 ) m ( (P - Qf f (P - R? (Q - R? ' (A - 6) 

where the notation is adopted from Using the finite temperature Feynman rules (see e.g. |35p in the 
Feynman gauge and taking advantage of the fact that the purely bosonic version ofT vanishes at O(e ) 0, 
we may write the diagrams of Fig. 1 in terms of the sum-integrals of Eqs. (|A.1|I - (|A.6|I . This computation 
is lengthy and requires the application of numerous tricks such as linear changes of integration momenta but 
is nevertheless of straightforward nature and is therefore not reproduced here. The result of the procedure, 
correct to O(e ), reads 

la = -(1- £ )^ 5 2 ^{X 1 (X 1 °-2X 1 )}, (A.7) 
/ 

h = 2(l-e)d A C A g i Y,{( I i-Zi)T + ^Mo,o}, (A.8) 

/ 

I c = (l-e)d A ^(7 F -i(7 A ^ 5 4 ^{4(l?-2^)f + (2 + e)^o,o-2eMo,o + 2M,-i}, (A.9) 

Id = -2(1- efd A C F g A Y J { (2? -2?) * % ~ 2I° 1 f + M ,o + M 1 ,- 1 }, (A.10) 

/ 

I e = -idAff 4 £{4(l + e)^[M/]^ fl ]Z2 + 2(l- e ^ (A.ll) 
fg L 

- 8 (ll\Mt]T'\ji g ] +X 1 1 K]r'[ M/ ]) - (1 -eJA/b.o^/,/^] - 2M,-i[M/,%] - 2Ar 2 ,_ 2 [ M/ ,/i ff ]|, 
If = -^d A C A g 4 Y / { 8I i^2-^if + M , -2M^ 2 }, (A.12) 
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I g = i^CVE{ 4 ( 6 - 5e ) X i° X i J 2 -(7-6e)I 1 T- (J -2e)AV -(5-4e).M_ 2 , 2 },(A.13) 
/ 

4 = -(3-2e)(l-e)d A C A g 4 ^{2X 1 X 1 X 2 -X 1 °f }, (A.14) 

/ 

where the sum-integrals Af n ,-n m Eq. (|A. 1 1|) depend on two independent chemical potentials /i/ and \x g 
through the respective fermionic momenta. The symmetry coefficients of the graphs have been taken into 
account here. 

Using the formula 

Zl = 1- UC V 47> X^ (A.15) 
9 3 (Air) 2 e 

for the renormalization coefficient of the gauge coupling, the unknown matching coefficient aE3 is now 
available in terms of the sum-integrals of Eqs. i|A.l|) - (|A.6(1 . One simply needs to add Eqs. (|A.7(1 - l|A.14|) 
together with the bosonic part of the strict perturbation expansion of the pressure, which can be found e.g. 
from Eq. (31) of We now turn to the actual evaluation of the unknown sum-integrals. 

Appendix B. Evaluation of the sum-integrals 

In this section all the sum-integrals encountered in appendix A will be computed. To do this we need to 
generalize the results of E] to finite /i, which includes determining the values of some new integrals of 
hyperbolic functions (see appendix C) as well as generalizing certain summation relations derived in [5]. 
The second but last section of this appendix deals exclusively with the problem of having different chemical 
potentials entering a single sum- integral, which is the case we encounter when evaluating the diagram I e . 
Until then it is assumed that all /i's inside a sum-integral are equal. 

B.l. One-loop cases 

The bosonic one-loop sum-integral I™ has been evaluated in 6 with the result 

jra = 2m -2n+l 7r m-2„+3/2 Tm -2„+4 ( A?Y ^ ~ ^ + ^ ((2n - m - 3 + 2e). (B.l) 

\7rT^ j r(?i) 

For the fermionic case we get after first performing a standard 3 — 2e -dimensional integral 

1 r(n-3/2 + e) a2 _ ~ ((2k + 1) vrT - i M ) m 



2"i-2n m-2n+3/2jrm-2n+4 f A 2 "\ T(n — 3/2 + e) 



i-3/2+e 



k ttT 2 J T(n) 

x U(2n - m - 3 + 2e, 1/2 - ip) + (-l) m ((2n - m - 3 + 2e, 1/2 + . ( B - 2 ) 
where the definition of the generalized zeta-function 

has been used. The results obtained for the relevant cases after performing an e-expansion are listed in Eqs. 

ESa - El. 
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B.2. Two- and three-loop cases 

The computation of the two- and three-loop sum-integrals in this paper is closely analogous to the /i = 
calculations of |H| and to the two-loop work of at finite fi. The general scheme is to first separate the 
diverging vacuum parts from the integrals and then to evaluate the rest in d — 3. This allows a straight- 
forward extraction of the singularities and also simplifies the determination of the finite parts considerably. 

The sum-integral T is not considered here, as its contribution will be observed to cancel between the 
different terms of Eq. (|A.11J| . 



B.2.1. Preliminaries 



Let us start by deriving some results for the bosonic and fcrmionic 'polarization' functions, 

1 



II(P) = 

MP) = 



qQ 2 {p + qY 
i 



\2 ' 



(B.4) 
(B.5) 



' {Q} Q 2 (P + QY 

which will be frequently used in the following computations. It is straightforward to verify that in the case 
of a bosonic external momentum P they can at d = 3 be written in the forms 



n(P) 
n f (p) 



T 



T 
(4^ 



dV 



(|po| +cothf)e- |po|r , 



where 



d 3 r 



Po 

f 



(B.6) 

,|r , (B.7) 

(B.8) 
(B.9) 

The latter formula, Eq. 1|B.7J| . is a generalization of Eq. (4.2) of [Hj and has been obtained using the 
Fourier-transform of the fermionic propagator, 



(M + cos(2/w)cschr) i 



Po 

2nT ' 
27rTr. 



d 3 q 



-(|2n+lj-irT-i/jsign(2n+l))r 



(B.10) 

(2tt) 3 q 2 + (( 2n + 1) ttT - ifi) 2 4ttT 

The large- P behavior of the polarization functions has been analyzed in [H] for the /i = case. A 
straightforward generalization of the computation to non-zero chemical potentials produces for Ilf 



n f (p) = n(°)(P) + n| T) (p) = n(°'(P) + n[ T u ) v (P) + o 



p6 



/P2 V 
\P 2 



Pi 



T 



d 3 7 



cos(2^r)cschr 



1 



1 



2/2 ; 
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360 



El. 
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2fL 4 

IT 



(B.ll) 

-|pok 



where e has been set to zero in the last 0^2^ -term. Here the coefficients read 



A) 



1 r(e)r 2 (l~e) /A 



(4tt 



.2— e 



r(2-2e) \T 



2c 



2 = 2 ^-3/ 2+ .r(3-2 £ ) /A 

H r(3/2-e) V? 1 



Li4-2 £ (-e 



Li 4 _2e(-e" 



(B.12) 
(B.13) 
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with Li denoting the polylogarithm function, and the relation 



LU(-e-) + LU(-e-)=-^-^*»-i^ 



(B.14) 



has been employed in deriving the last term of Eq. (|lj.llfl . The first term of Eq. (|lj.llf> . denoted by II^ ^, 
is the vacuum (T = fi = 0) part of the polarization function, and is the same in the bosonic and fermionic 
cases. Following the notation of the vacuum subtracted, or matter, part of the function is denoted by 

(T) 

LTj , even though here it should not be confused with the finite temperature part of Ilf obtained simply by 
subtracting the T = piece from the polarization function. 

B.2.2. T 

In terms of the function Ilf and its large- P expansion, the sum-integral T can clearly be written as 



r 



i 

T 
(4^ 



MP) 



d 3 r- 



1 



P P 2 
cos(2/if)cschf 



p2 



22j 



1 

P 2 



2p,' 



I Poh 



0(e) 



The first two terms are trivial to evaluate using Eq. (|B.1|I , whereas the last one gives at d = 3 



(B.15) 
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(Air) 2 % P 2 



T 2 
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(4nT 
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cos(2//r)cschr — \ 
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-bol' 



cos(2/ir)cschr - 



2/2* 



dr 



cothr 



cos(2/ir)cschr 



2/2' 



(B.16) 



Here the final integral is UV convergent but has an IR divergence due to the zero mode of the bosonic 
frequency sum. Using the results of appendix C it can, however, be evaluated analytically and produces a 
finite result. This is due to the fact that the divergence in the r-integral comes from a term of the form 
J drr Q , which vanishes under dimensional regularization. Adding up the different pieces and expanding 
in e, we get the result, Eq. (|B.73(1 . for the whole sum-integral. 



B.2.3. W ,o and M ,o 

The application of Eq. 1B.11|) separates the evaluation of No,o into three pieces 



M),o 



^ (nfV)) 2 = ^ ((n< )(P)) 2 + (nf ) (p)) 2 + 2n(°)(p)nf ) (p)j , (b.17) 



of which the first one is again trivial. The second produces in analogy with the above two-loop calculation 
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(B.18) 
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where the resulting one-dimensional integral is analytically calculable. 
The third term of Ao,o can be divided into three parts 



+ -^r f p ( n f T) ( p ) ~ ng v (P)) + |f n(°) (p) nf^P) 



(B.19) 



(B.20) 



Here -K2 and K3 are straightforward to obtain using the previous results, but the first term requires careful 
consideration. Taking into account that 



^ (nf^-n^p)) =o( £ 



and using the result of (Hj, 



{2nf 



4vrA 2 
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(B.23) 



V360 3 3 

The remaining integral is again of analytically solvable type, and after adding together all the different 
pieces we get as the result Eq. (IB.74|I . 

The computation of A'Jo.o proceeds in a similar fashion. Eq. (|B.17|I is replaced by 



M j 



n (0) (p)j +n (T) (P)nf \p) + n (0) (P) (n (T) (p) + n f (T) (p)j , (B.24) 



where the first term is trivial and the last two cross-terms are available through Eq. (|B.20(1 above and Eq. 
(D20) of UJ. For the remaining piece one obtains 



n<- T )(p)ny>(p) = 
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whose straightforward evaluation leads to the final result, Eq. (|B.75|) 
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B.2.4. Wi, i and Mi-i 

It is easy to see that the sum-integrals and A^i.-i can be written in the form 

Q-R 



= 2X 1 °T - 2 

= 2X 1 °T - 2 
Following 6 and defining further 

M P ) 

MP) 



we quickly verify that H 3 and H3 read 
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f p -ps MP) MP) 



h 3 
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I Z 3 +4^0,0, 



1 ^-5?) 9 2j-i 



(B.26) 
(B.27) 

(B.28) 
(B.29) 
(B.30) 
(B.31) 

(B.32) 
(B.33) 



This leaves only the simpler sum-integrals I 3 and I3 to be evaluated. 

Due to Lorentz invariance and the orthogonality of J'^{P) to bosonic P and J^(P) to fermionic P it is 
evident that these functions can be written in the forms 



MP) = 



MP) 



p2 

2 

F 2 



Po 

p2 



Pjio(P). 



which, when plugged into Eqs. (|B.30fl and i|B.31|) . produce 

A crucial simplification in the evaluation of these sum-integrals occurs, as one notices that they both are 
actually finite. This is due to the fact that at large P J^(P) and J^{P) behave like O (l/P 2 ), as can be 
straightforwardly verified. We may therefore set e = in the expressions for f and jo, which eventually 
yields 



(B.34) 
(B.35) 

(B.36) 
(B.37) 
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sign(Repo)- (B.39) 
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After substituting Eqs. l|R38)l and l|B~39|) to Eqs. (|R36)l and l|B~37) . doing the Po-sums and performing 
the remaining p-integrals using the relation 



d 3 p 1 sin pr sin ps 



(2tt) 3 P 2 pr ps 



we finally get 



i; = - 
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(B.40) 
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+ i((l + 4/2 2 )r + 4 i/ 2)) 1. (B.42) 



1287T 2 

These integrals are clearly both UV- and IR-finite and are readily evaluated. The final results for the 
sum-integrals, obtained after adding up all the different pieces, are given by Eqs. <|B.76() and 1|B.77(1 . 

B.2.5. W 2 ,-2 and M-2,2 

Let us define, once again in analogy with [S], the modified gluon polarization tensors, 

(2Q + P) M (2Q + P)„ 



LV(P) = 25^1" 



Q 2 (Q + P) 2 
(2Q + P) M (2Q + P\ 



'{q\ Q 2 (Q + py 

Using these we further define two new sum-integrals by 



sqed 



T IU = % -^An^(F)An f ,^(p) 
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(B.43) 
(B.44) 

(B.45) 
(B.46) 



where A/ (P) = / (P) — / (0)5 po . It is then straightforward to verify that determining Ni^-i and M.-2.1 
can be reduced to computing a set of simpler sum-integrals 



where 
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\ ill* + \ + \ M ^ - J " T + ( 2 + 2e ) 



The evaluation of PjZ d and PZ ed is fairly easy, since an elementary calculation verifies the validity of Eqs. 



d 3 - 2t p 1 - 
(2tt) P 

d 3 ~ 2e p 1 
(2nf- 2e 7 



S M „ (o) n f;M „ (po = 0, P ) + n f ,^ (o) (p = o, P ) 



(B.47) 
(B.48) 

(B.49) 
(B.50) 



(F6) - (F10) of [Hj also when n f 0. This gives at order O (e) 



n[ T) (p)) +4(d-2)(2y) x 2 u , 



(B.51) 
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|£ ^AflP?(p)Ang,(p) 
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^ p (n»>(P)) 2 , 



(B.55) 



where the superscript (0) again signifies the vacuum part of the function in question and f^ T > = f — f( ' 

The relations (|fj.51|l - (|B.55|) clearly reduce I^ ed and I b sqed to well-known sum-integrals and leave only 
^Iqed an d visaed to ^ e evaluated. Noticing that for P = (0,p) 



4P„ 



Q 2 -<& 



3 - 2e r {Q} Q 4 
^-(( 2 -4^ + 4T|) 



(B.56) 



and similarly 

p^n M „(o) = fyif^o.p) = p M n^(o, P ) = o, (B.57) 

we may decompose the tensors into their transverse and longitudinal parts (see e.g. 123). This enables us 
to write Eqs. i|B.49l) and (|B.50(I in the form 
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and the use of the identity I^l^-i = ^f/^i = — 1/2 + e now leads to 
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In Eqs. 1B.60|) and 1B.61|) the p- integrals can be performed by introducing Feynman parameters x and y 
to combine the different factors in the denominators. A straightforward calculation produces 
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where the function f is defined by 
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The use of this result in Eqs. (|B.62|I and l|B.63l) now gives at order O (e) 
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and together with the other parts of A/2,-2 an d this yields Eqs. I|B.78|) and l|B.79|) as the results. 



B.2.6. The diagram e 

Eq. (|A.11|) . the contribution of the diagram e to the strict perturbation expansion of the QCD pressure, 
contains sum-integrals, which are functions of two independent chemical potentials. Due to cancellations 
between the different terms of I e it is convenient not to deal with all of them separately, but to treat the 
whole diagram as a single entity. 



Using trivial generalizations of Eqs. 1B.47|) and l)B.51|l we get for the diagram 
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where every term is of an already familiar form with the exception of the one containing the r-integral. This 
integral can, however, also be straightforwardly evaluated using the results of appendix C, which eventually 
gives us Eq. I|B.80|) as the final outcome of the graph. 



B.3. The results 

The final results for the sum-integrals introduced in appendix A and evaluated above read 
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Appendix C. Evaluation of the hyperbolic integrals 

As a consequence of keeping [i finite, many of the one-dimensional integrals encountered in this paper differ 
from the ones of |S] . In addition to the ordinary hyperbolic cases we need to evaluate integrals of the type 
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where n and p are non-negative integers and z and f3 real numbers. This is accomplished by repeatedly 
applying the relations 

coth 2 x - csch 2 x = 1, (C.2) 
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which can be straightforwardly derived. As in [jy, UV-divergences in the individual terms of converging 
integrals are regulated by introducing a factor x s in the integrand and in the end taking the limit 6 — > 0+. 



Appendix D. Properties of the functions N 

In section 2 the functions H were defined by the formulas 

C(x,y) = d x ((x,y), (D.l) 

X(n,w) = C(-n,w) + (-l) n+1 ('(-n,w*), (D.2) 
N(u>) = $(w) + *(«;*). (D.3) 

In order to analyze the behavior of the sum-integrals of appendix B at different values of fi and T, we 
need to expand these functions in the limits of small and large /2. The results of such expansions, obtained 
straightforwardly using among other things the integral representations of the zeta and digamma functions, 
read 
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Appendix E. The diagram and the plasmon term at T = 

In section 4, where the zero-temperature pressure of QCD was computed to 0(g 4 ), we only quoted the 
results for the diagram I' e and the plasmon term. The corresponding calculations will be performed in this 
appendix. 

E.l. The diagram I' e 

The one-loop gluon polarization tensor depicted in Fig. 4.b. has in the Feynman gauge the expression 

(ll^rV) = —2g 2 Tp5 ab (ffiV + (P,P V - P\ v ) U f (P) - g£ {2Q ~J^ Q p ~ P)v ) ■ (E.l) 

Denoting its vacuum (T = \i = 0) part here by (n^) ab (P) | vac we obtain after a straightforward compu- 
tation 

(n f r 6 (p)| vac = ^^(^'(p^-pV). (E-2) 



where the coefficient ^4 can be shown to have the the e expansion 
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Using this expression we easily obtain at T = 
I' e = A(d-2)g 4 Ad A T^A 2e 
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and only the integral T remains to be evaluated. 

The most straightforward way to tackle the computation of the new integral is to proceed as in the case of 
the sum- integral 7~, while from the beginning on neglecting terms that vanish at T = 0. Using Eq. (|B . 1 1|> 
we get 
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where the remaining P-integrals can be evaluated using previous results. The calculation of the r-integral 
then leads after setting the temperature to zero to 
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and plugging this expression to Eq. (|E.4|) gives Eq. 15. 2|) as the result. 



E.2. The plasmon term 

The evaluation of the plasmon contribution to the T — pressure corresponds to summing over all the ring 
diagrams of Fig. 4.c. starting at three-loop order. These graphs contain as loop insertions the fermionic 
part of the vacuum subtracted one-loop gluon polarization tensor, which we denote by (AII^) . The 
summation was originally performed in |2] and will be reproduced here following to a large extent the 
treatment of I38| . The computation begins from the derivation of an integral representation for the 
polarization tensor and then proceeds to the evaluation of the actual plasmon sum. 

Using Eqs. (|E.1(1 and (|E.2|) it is easily seen that the vacuum subtracted polarization tensor is orthogonal 
to the external momentum, i.e. 
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Setting e = and applying this relation together with rotational invariance we may then write the tensor 
in the form 
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where P stands for four-vectors and p for three-vectors as before and the polarization tensor has been 
divided into two parts proportional to orthonormal projection operators. The components of the tensor 
appearing explicitly in the result are available through an application of the Residue theorem, which after 
a straightforward calculation leads to the results (see appendix A of part II in [H]) 
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and performing the corresponding integrations in Eqs. ijE.lOf) and i|E.ll|) . we now obtain for the polarization 
tensor 
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where P this time stands for the norm of the corresponding four-vector and where flavor sums have been 
suppressed. The last forms obtained here are, however, inconvenient to work with. It is on the other hand 
easy to confirm by a direct integration that the simple integral representations 
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produce exactly Eqs. (|E . 13|> and (|E. 14(1 . which suggests the use of the Sommerfeld- Watson integral formula 
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When applied to Eqs. i|E.15(l and 1E.16|) . the Sommerfeld- Watson formula gives 
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where the x and y integrals can now be factorized and performed separately. The result first derived in 
is that one has obtained a compact integral representation for the polarization tensor 
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where flavor-sums have been re-introduced and the functions r„ are defined by 
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We are now ready to evaluate the sum of the diagrams of Fig. 4.c, which we begin by first defining (again 
in analogy with 0) 
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Using the orthonormality of the projection operators appearing in Eq. 1E.9() and keeping in mind the 
symmetry factors for the ring diagrams, this straightforwardly leads to the expression 
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which we now must examine. 

The integrals over the three-dimensional angular variables are trivially performed. This gives for the 
infrared sensitive part of the plasmon contribution 
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where the P 2 -integral has been evaluated in the last stage. From Eqs. (|E.13(I and 1E.14|) we furthermore 
have 
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Performing the remaining 0-integral produces now 
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with the constant <5 defined in Eq. (|5.5|l of section 4. 
The second part of the plasmon term is defined by 
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To order g 4 it can be evaluated by simply expanding the logarithms in Eq. I|E.26|) in powers of g 2 , since 
there obviously will be no g 4 ln g contributions originating from its expression 



P\ 



dP 2 P 2 



tt/2 



isin 2 6 s In 



P 2 +A 1 (P) 



+ 2 In 



P 2 +A 2 (P) 
P 2 + A 2 (O,0) 



W Jo Jo I L^ 2 + Ai(0 

^ (Ai(P) + 2A 2 (P) - A x (0, 6) ~ 2A 2 (0, 0)) - ^ p2 j 4/A 2 ( A ?(°> + A 2(0, 

t/2 



- 2<iFr dp2 r^ sm2 Hi( A ' (p,+2A = (p) 



1 



1 



A 2 (O,0) + 2A 2 (O,, 



v P 2 + 4/x 2 P 2 

Plugging in the contour-integral representations of the A-functions 
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and performing the P 2 -integral with an IR cutoff Ai] /x 2 , we now obtain 



2(2tt) : 



4?7/X 2 



f 

P 2 



fg 

lo °- e dz_ 
2nd 

tt/2 



tt/2 



>sm 







f 



n 2 V^ 2 + 4/x 2 P 2 



-i- ) (A 2 (0» + 2A^(0, 



(E.38) 



— lOO — € 



zo °- e dz^_ 
, 2ttI 



— IOC — 6 



(A{(z,0)A?(z',0) + 2A|(z, 0)A| (z', 0) 
2(2^)3 E/ # S in 2 0|-^(A 2 (O^) + 2A 2 (O,0))lnr / 

. ^ 7 (A{(z,0)A?(z',0) + 2A^(z,0)Af(z', 



P 2 
44 



(E.39) 



fg 

lco ~ e dz 



— 'ioo — e 



2m 



4 °°- £ ' dz' 7f+ z ' 



m2 



The ^-integral in the first term is easy to perform, but is left in its present form for now. Instead, we 
deform the integration contour of the z' integral in order to be able to proceed to the limit 77 — > 0. Denoting 
Residues by 5ft and encircling the poles at z' = and z' — —z to obtain a contour-integral, which vanishes 
in this limit, we get 
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The z-integral of the first term is quickly evaluated by closing the contour on the right half-plane, where 
the only pole is located at the origin z = 0. The substitution rf = 1 + zlnr/ + 0(z 2 ) then produces 
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where we have used the relation 
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a trivial consequence of the definitions of Eqs. (|E.36|I and (|E.37|I . 

The evaluation of the plasmon term has now been reduced to the computation of the two integrals pg 
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the double residue appearing in p 3 ' 2 is rapidly evaluated, and the computation of the subsequent ^-integral 
produces 
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For p 3 ' the (^-integration can, on the other hand, be immediately performed. This gives 
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where the remaining contour integral is straightforward to evaluate by applying the Residue theorem. 
Assuming first /i 2 > y? g we may close the integration contour by a semi-circle on the right half-plane, where 



the integrand obviously has poles located at z — k, Vfc <E N. Denoting y? q j \A = (3 we get 
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where the infinite sum has been evaluated using standard formulae. For fij < fi 2 <^> (3 > 1 the computation 
proceeds in an analogous fashion the only difference being that the integration contour is now closed on the 
left half-plane. The result of this calculation is simply 
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and the application of Eqs. I|E.49|) and (|E.50() in Eq. IE. 481) gives now 
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Adding Eqs. (|E.33(1 . I|E.48|) and (|E . 5 1|) together we finally obtain Eq. I|5.3|l as the result for the whole 
plasmon sum. 
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